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Abstract 



We develop a formulation of the Baade-Wesselink method which uses the Fourier coefficients 
of the observables. We derive an explicit, analytic expression to determine the mean radius 
from each Fourier order. The simplicity of this method allows us to derive the uncertainty 
\q in the mean radius due to measurement errors. 

Using simulations and a recent dataset we demonstrate that the precision of the radius 
measurement with optical magnitudes is in most cases limited by the accuracy of the mea- 
surement of the phase difference between the light and the color index curve. In this case it 
^ is advantageous to determine the inverse radius, because it has normal errors. 

y—i 1 Introduction 

The Baade-Wesselink (BW) method of measuring radii and distances of Cepheid variable 
h stars (Baade 1926, Wesselink 1946) is potentially very useful to establish the zero point of the 
Cepheid period-luminosity (PL) relation. Gautschy (1987) gives a good review of the status 
of the BW method. Gieren, Barnes & Moffett (1993) have used a sample of 100 Cepheids to 
derive the PL relation, which agrees within the uncertainties with the ZAMS-fitting distance 
Q scale. Laney & Stobie (1995) have used infrared photometry to derive a period-radius (PR) 
q relation that is steeper than previous results. Pipepi et al. (1996) use an improved two color 
£j index method to find a much shallower PR relation. All the results mentioned above are 
& based on different BW methods, photometric bandpasses and data sets. Laney & Stobie 
(1995) showed that the choice of photometric bandpass as well as the BW method affects 
the resulting PR relation. 

This is the first paper of a series to improve the PR relation. In this paper we introduce our 
BW method and discuss the sources of measurement errors in BW radii. In subsequent papers 
we will apply our BW method to Cepheids and RR Lyr stars in the Galaxy and the LMC. 
Estimating the uncertainty in the radius of each individual Cepheid due to measurement 
errors is a problem in all realizations of the BW method. Both photometry and radial 
velocity measurements have errors that lead to an uncertainty in the BW radius. Knowing 
this uncertainty is important to be able to look for systematic errors. Good error estimates 
are also essential for determining the analytical form of the PR and the PL relation. Without 
an error estimate for each individual Cepheid radius or distance the quality of a PL or PR 
fit cannot be determined. Even though linear PR and PL relations are most widely used, 
there is still uncertainty whether or not this is the correct analytic form. Burki and Meylan 
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(1986) find evidence for a break in the PR relation at a period of 10 days and Parsons (1972) 
fits a quadratic PR relation. 

Balona (1977) includes the estimated measurement errors to derive maximum likelihood 
radii. The uncertainty Balona (1977) quotes for each stellar mean radius is the standard 
deviation of the individual radius determinations at each available photometric point. These 
errors include uncertainties in the phasing of photometry and radial velocities as well as 
possible problems with the assumption of a single valued color-surface brightness relation. 
The effects of the measurement errors cannot be separated from these other errors. Another 
problem with propagating measurement errors in the maximum likelihood method is the 
required interpolation of the velocity curve. This makes it difficult to assign velocity errors 
to the phase points used in the solution for the radius. 

It is advantageous to use phase-integrated quantities to measure the BW radius because 
no interpolation is required and integrated quantities are less sensitive to noise. Caccin et 
al. (1981) introduced a BW method that determines the radius from an implicit, phase- 
integrated equation. In order to get phase integrated quantities they fit Fourier expansions to 
the data, which can easily be integrated. To estimate the influence of the measurement errors 
on the radius Caccin et al. (1981) did Monte Carlo simulations of several noise amplitudes. 
Using a certain £ Gem data set they estimated the radius to be uncertain by about 5 % . 

Another phase-integrated BW method uses the Fourier coefficients of light, color and velocity 
curves (Balona and Stobie 1979). Assuming a certain relation between the shape of the 
velocity curve and the shape of the color curve, the Fourier components of the light, color 
and velocity curves are used to determine the mean radius. The scatter in the resulting 
PR relation for 60 Cepheids using this Fourier method is larger than the scatter in the PR 
relation using a maximum likelihood method (Balona 1977). Thus the assumptions of this 
Fourier method might not be a good approximation for all Cepheids. 

In this paper we present a new BW method using the Fourier expansions of the measured 
quantities to determine the radius. An advantage of this method is that measurement errors 
can easily be propagated. We go on and use this to discuss the influence of errors on the 
estimated radius. 

2 The Fourier Baade— Wesselink Method 

The radial pulsation of Cepheids provides a direct means of determining the absolute radius 
of each Cepheid via any of the existing variants of the classical BW method. The method 
is based on observing the radial motion of an emitting surface through the Doppler shift 
in its spectrum, as well as its flux and temperature variation. The radial velocities, when 
integrated appropriately over a certain time (or cycle), yield the linear displacements of the 
moving stellar surface. If the temperature were constant, the flux variation would reflect 
the change in radius of the emitting surface, and yield radius ratios at the same times. 
Combined, these two radius variations would lead to a solution for the mean radius. 
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However, as it pulsates a Cepheid changes its surface temperature. Both radius and tem- 
perature changes affect the total luminosity, which we observe as the Cepheid light curve. 
The effective temperature at the surface of any Cepheid is between 5000 and 7000 K; the 
variation due to the pulsation is usually about 1000 K. The emergent spectral distribution 
for that temperature range peaks in the optical region. As a result, the optical light varia- 
tion is dominated by the surface brightness variation, while we need to determine the small 
variation caused by the change in radius. A common way to separate the surface bright- 
ness variation is to use the observed variation of the photometric color index as a proxy for 
the surface brightness. Then we have three observables: two light curves at two different 
wavelengths (bandpasses), and a velocity curve; we solve them simultaneously for the mean 
radius. 

In this section we derive an equation that gives the mean radius of a pulsating star in terms 
of the Fourier coefficients of these three observables. This equation is very useful because 
it allows easy propagation of the measurement errors and gives insight into the main source 
of error in the radius. We assume that we have data with good coverage over phase, <fi, of 
the observables and that the period of the Cepheid is known. The starting point of our BW 
solution is the relation between magnitude, surface brightness and radius 

/ R \ 

M = S - 51og( — ) + const, (2.1) 

\Rq/ 

We assume that the surface brightness can be written as a power series of the tempera- 
ture observable, T(<p), S(cf)) = Co + J2j=i c i^' 7 (^)- The temperature observable might be a 
color index, or a temperature sensitive line ratio. In this paper we define the temperature 
observable to be the B V color index. 

Then all observables, which are measured radial velocity V, magnitude M and color index 
T, can be Fourier expanded to give 

m 

X = X + ^X fc cos(27rfe((/) + (/)f)) (2.2) 
k=i 

where <fi is the phase ranging from to 1. 

To write eq. (2.1) in terms of the Fourier coefficients of the observables we need to express the 
radius in terms of the velocity. We allow for a phase-dependent conversion factor (p-factor) 
from radial to pulsational velocity, p((f)), thus 

R(4>) = Ro + Ar(0) = R + P t p^') (v (</>')- 7) d<j>' (2.3) 

Here Rq is the unknown mean radius, P the period of the Cepheid and 7 the systemic 
velocity. We assume that p(4>) can also be Fourier expanded, so that we can write 

Jo p(<j>)V(<j>)d<j> 1 v-^t/ to -(Ay Aj>\\ (OA\ 

Jo pwW Zpo 3=1 
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Even though we set up the equations very generally, we will now solve them for the simplest 
case. This demonstrates the concept and keeps the algebra simple. Later we will discuss a 
more general solution. In the simplest case the p-factor is constant over phase and the color 
index-surface brightness relation is linear. Thus p = po, 7 = Vq, S = cq + c\T and 

k=i 

Substituting all Fourier expansions, equation (2.5) and the linear color index-surface bright- 
ness relation into eq. (2.1) we have the following equation. 



M + J2 M k cos(2ttA;(^ + </>f )) = c + ci(T + T k cos(2tt^ + $))) 

k=l k=l 

-51og(l + f;|^sin(27rM^ + ^))) - 5 log (-^) + const. (2.6) 

The key step in our BW formalism is to separate out the Fourier orders from each other. 
This we can do if we linearize the logarithmic term in eq. (2.6). We use ln(l + x) ~ x, if 
x <C 1. In our case x = Ar/Ro, which is typically less than 0.2 for PopI Cepheids, thus 
the approximation might be acceptable for these stars. We will check the validity of this 
approximation below. 

After approximating log(l + x) by xlog(e) we use the orthogonality of sin and cos to break 
eq. (2.6) up into a system of equations. We multiply eq. (2.6) by sin(27rj</>), or cos(27rj</>) 
and integrate over phase to find 

„ 5 \og(e)poPVj cos(2n jd>Y) 
Mj sin(27rj0f ) = Tj sin(2-. / o / / )r l / + ^ h 



2irjR 0j 

„ 5 \og(e)poPVj sm(2nj(f)Y) 
Mj cos(27rj^ M ) = Tj cos(27rj^ T ) Cli 2njR j (2 ' ?) 

Note that these equations are not a series of approximations. Each Fourier order contains 
all the information, thus in the absence of noise i?oj = Roi an d c\j = c\{. Combining the 
two equations, the mean radius, Roj and the slope of the color-surface brightness relation, 
cij, are given by 

R 0j = 0.3456 j~ pqP J . ) n ' M J " 



M j sin(2 7 rj(^f - <f>J)) 



M i cos(27rj(^-^f)) 
r iC os(27rj(^-^)) 



In the absence of noise, each Fourier order j will give the same values of Roj and c\j. For 
real Cepheid light, color and velocity curves the values of Raj and c\j will be affected by 
errors and thus be different for each Fourier orderx. 
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2.1 Effects of Linearization 



The correct radius term of eq. (2.1) is given by 



h{<j>) = 5 log ( 



Rq 



(2.9) 



This must be compared to the linearized radius term, given by 



/ 2 (<£) = 51og(e) 




(2.10) 



We found that the Fourier amplitudes of typical Cepheid radius curves are reduced by the 
linearization (see below). This decrease is smallest for the first order and increasingly larger 
for the higher orders. The Fourier phases do not change at all between the linearized, fi 
and the full radius term, f\. 

The reduction in the amplitude of the radius curve leads to a smaller mean radius. This 
effect increases with large velocity amplitudes, the steepness of the piston phase and long 
periods (large Ar). The effect of the linearization on the determined radius can be removed 
by a correction factor given by the ratio of the Fourier amplitudes of the fi,2(4>) curves. 
The procedure for determining this correction factor is the following. We first solve eqs. 
(2.8) to get the mean radius, Rq. With this mean radius we determine and their respec- 
tive Fourier coefficients. The ratio of the ith order Fourier amplitudes of f\ and fi is the 
correction factor for the mean radius determined from the zth Fourier order. 

Below we will examine the effectiveness of this procedure with simulations. 
2.2 Measurement Errors 

The simplicity of eq. (2.8) allows direct propagation of the measurement errors to estimate 
the uncertainty in the derived radius. 

2.2.1 Errors in the Fourier Fit 

For each of the observables X, where X = V, M or T, the Fourier amplitudes, X n , and 
Fourier phases, , as well as their respective uncertainties are derived from a least squares 
fit of the velocity, magnitude and color index curves to Fourier expansions of the form 



We assume that the period of the Cepheid is well known and that the velocity, magnitude 
and color index curves have been phased whith this period. The coefficients X^ and (j)^ of 
eq. (2.2) are related to a\. and by 



m 




(2.11) 



k=l 




(2.12) 
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and 

arctan 

iX _ 



V a k ) 



k 2irk 



(2.13) 



The variances and covariances of and (p^ are given by 

0'k a l cr a k ai + ^kh°b k bi + a kh^a k bi + ^l^k r a t b k 



°X k X l 



4**k 1(4 + + 
o-k^akbi - bkbiab kai - a-khvakak + bkdia bkbl 

2M + bf) x /4+~q 



= : . r= r I 2 - 14 ) 



The uncertainties, c a & a & are given by the covariance matrix of the least squares fit. 
2.2.2 Errors in the Radius of each Fourier Order 

The radius is determined from the coefficients Xf. and <jr£ by eq. (2.8). The errors in the 
radius are given by X^ and and their variances and covariances. We assume that there 
are no correlations between the Fourier coefficients of velocity, magnitude and color index. 
The variances and covariances of the mean radius are given by 

CRq-Rq- Oy.y. OMM 

Rq * Rq J . = ~vf\/~ + M-M- + ^^^^Y^J + frm(i)fr m (J)°<l>M<l>M + fr t (*) fr t C?>#>T + 

0') _ 



with 



3 v% J lv± 3 



f rv (j) = -2njtan{2nj{<l>Y - $)) 



/ rm (j) = -2 7 rjcot(2 7 rj( ( /»f -<#)) 
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COs(27Tj(^-^ M )) 

}r ' 0) " '^sin^^f - p cos(27rj(^ - p ( ~ * L6) 



Likewise the variances and covariances in the slope of the color-surface brightness relation, 
ci, are given by 

&ci ci ■ &T T- &M M- 

— ~ = + m T7 + fcvi^fcvU^Yd^ + kWkO>^^ + fc t {i)fc t {j)(J <t> T < hT+ 

C\iC\j -Li-Li Milvlj v * v 3 Y i r j r i 

/rt(«)_ , /r t (j) k(') /r m Q") , 7 n 



with 



fc v (j) 



sin(27rj(<ftf - ffi) 



fc m U ) 

fctti) 




(2.18) 



2.2.3 Do the Errors Follow a Normal Distribution ? 



The radius, as given by eq. (2.8), is the ratio of two terms, Roj oc Xj/yj. 



Xj = y j cos(27rj((/) i - )) 
Vj = Mj 8 m{2n j{<f>f -<f>J)) 



(2.19) 



Xj and its error is mostly determined by the velocity curve and yj and its error by the 
photometry. If the probability distribution functions of Xj and yj are known, Pji(xj) and 
Pjziyj), then the probability distribution of i?oj is given by 



Below we will use simulations to demonstrate that pj\ and pj2 are well approximated by 
normal distributions. In this case the above integral can be written in terms of the error 
function, but still must be evaluated numerically. 

To get some insight into the probability distribution of Roj we consider two limiting cases. If 
the error in Xj is dominating the error in Roj, the probability distribution of i?oj is normal. 
Eq. (2.15) then gives a correct estimate of the width of the normal probability distribution. 
If the error in y is dominating the error in Roj, the probability distribution of i?oj is n °t 
normal. But the inverse Radius, R^j follows a normal distribution. In this case the error 
given in eq. (2.15) is not a good description of the true error. In this case we can instead use 
the inverse radius, Roj, because its errors are normally distributed. The error of the inverse 
radius is easily calculated from eq. (2.15). 

In Figure 1 we show the probablity distributions of the radius for a normally distributed 
inverse radius with a a mean of 0.01 and a a of 0.004. The probability distribution of the 
radius is asymmetric and the errors are different for positive and negative deviations. 

2.2.4 Combining the Radii 

The best mean radius is obtained by minimizing the likelihood function 




(2.20) 



L = l[P J (R 0j -R ) 



(2.21) 
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Figure 1: A Normal distribution of R 1 leads to an asymmetric distribution in R. The distribution of 
radii peaks before the true radius and then there is a long tail of large radii. 



In the case of normally distributed radius errors the final radius is calculated using 



Ro 



Ej °rI,R, 



0j j 



The error is given by 



(2.22) 



(2.23) 



The same equations can be used for the inverse radius, if the probability distribution of the 
inverse radius is normal. 

2.3 More General Solutions 

In this chapter we deal with all the terms we omitted when deriving eq. (2.8). 
2.3.1 Phase dependent p factor 

We now assume that we have the p-factor as a function of phase. The p-factor can be derived 
using models of Cepheid atmospheres (Sabbey et al. 1995) or using the measured widths 
and asymmetries of the spectral lines. With a phase dependent p-factor eq. (2.5) becomes 



V 



Ar = P ( po ^ sin(27rj(^ + <%)) + (V - 7 ) £ s in(27rfc(</> + $)) 

j=l 713 k=\ * 



, S m(27r((j-k)<f> + j<f>Y -k<%)) S m(27r((j + k)<f, + j<f>j + k^ k )) ^ 
l^ V i Pk \ 4n(j-k) + 4n(j + k) ) 



j,k—l 



The radius becomes in 1st order 
3456 P t 

R ° = T« ■ (o (AM MTW ( PoVl C0S (M<^ " <t>D) + (V - 7 )pi cos(2tt(^ - </>f )) + 



£ ^ cos( 27 r(0- + 1)^+1 " ^ " #)) + 
j'=i 

£ ^ cos( 27 r((, + - ^ - #))) (2.25) 



2.3.2 Non— linear Color Index Surface Brightness Relation 

Allowing for non-linear terms in the color index-surface brightness relation eq. (2.7) becomes 



Mj 8in{2n j<f>f ) 
Mj cos(2n j<j)f ) 

where and </>J^ ^ are the jth order Fourier components of T k . In order to solve this 
system of equations we need to use higher Fourier orders. This requires high quality data, 
because if the higher orders contain little information, the solution will be dominated by 
noise. Using the first and second order Fourier components will allow to include terms up 
to T 3 . This should be enough for most practical purposes. 

3 Using Simulations to Asses the Effects of Errors 

In order to demonstrate the effects of linearization and measurement errors on the radius 
we studied simulated data sets. The magnitude and velocity curves used for the simulations 
were derived from Fourier fits to the velocity and the magnitude curve of U Sgr, obtained 
by Bersier et al. (1994). For each phase point the Fourier fit allows us to "measure" the 
magnitude or the velocity. In order to exclude all systematic errors we assumed a mean 
radius and constructed a matching color index curve that solves eq. (2.1). We assumed 
either 20 or 40 points in both velocity and photometry. The phases of the data points were 



£7f sin( 27 r^fV + 



k=l 



cosine* 



k=l 



5 \og(e)pPVj cos(27r j(f) - 

27TjR 

5 \og{e)pPVj sin(27rj0p 

27TjR 



(2.26) 
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chosen randomly, but we required that there be no large gaps. The number of Fourier orders 
fitted for was determined from the largest phase gap, A max , in the data to ensure that the 
Fourier fit does not blow up between two data points. We chose the number of Fourier orders 
to be the smallest integer larger than (A"^ — l)/2. The typical number of Fourier orders 
fitted for was 3 - 5 for 20 data points and 6 - 8 for 40 data points. 

An important quantity in determining BW radii is the phase difference, Acfij, between the 
jth orc [ er Fourier phase of the color index, T, and magnitude, M, and its error, o^.. This 
phase difference would be zero if the radius were constant. As the change in radius is small, 
Acj)j also is small, because most of the change in brightness is mostly due to the temperature 
change for an optical bandpass. If infrared (IR) magnitudes and colors are used A<pj becomes 
larger. 

Because A<pj is very small and in the denominator of eq. (2.8), the error in the radius is 
dominated in most cases by g^.. A small change in Acfij leads to a large change in radius. 
Thus the inverse radius has normal errors, whereas the probablity distribution of the radius 
itself is asymmetric (Figure 1). Therefore the average inverse radius of 1000 simulations with 
typical photometric and velocity errors gives a very good estimate of the true inverse radius, 
but the average radius is significantly larger than the true value due to a few outliers. 

The insights that we gained into the errors of BW radii are not restricted to our Fourier 
method. To demonstrate this we used a multidimentional minimization routine to solve eq. 
(2.6) simultaneously for all orders. The result of this procedure is a radius, derived from all 
Fourier orders, with no linearization. The radii derived using this method are very similar 
to the results of the Fourier method, but there is no information about the size and the 
behaviour of the errors. 

In Figure 2 we show the distribution of inferred radii from 1000 simulated Cepheid data sets 
with typical measurement errors. We used a straight multidimentional minimization routine 
to derive the radii. The distribution of radii is asymmetric. Because the errors are dominated 
by Acfij, the errors in 1/Rq are normal. The distribution functions are very similar to those 
in Figure 1. 

3.1 Biases due to Measurement Errors 

We used measurement errors with normal errors for which we give the width, a, relative to 
the peak to peak amplitude of the Cepheid. The velocity errors we used have a a of 0, 5, 10 
and 20 % of the amplitude of the velocity curve. For the photometric errors we used a a of 
0, 1, 3 and 5 % of the amplitude of the light curve. We assumed that the error in the color 
index is 0.75 of the photometric error. 

We simulated data of two different periods, 2 and 20 days, for which assumed a mean radius 
of 28.8i? and 97.6i? respectively. For a period of 2 days and our choice of velocity and 
magnitude curves A<pi = 0.029, and for 20 days A<pi = 0.073. A<pi becomes larger for the 
longer period, because At/Rq becomes larger. The radius corrections due to linearization 
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Figure 2: Distribution of the inferred radii from 1000 simulated Cepheid light, color and velocity curves. 
We used a multidimentional minimization routine, photometric errors of 3 % and velocity errors of 5 %. 
The inverse radii are well fitted by a Gaussian, indicating that the photometric errors are dominant. This is 
confirmed by applying the Fourier method to calculate the errors. The straight average of the radii is 10 % 
larger than the true radius. 



are small; the uncorrected first order radii are smaller by 0.4 % for a period of 2 days and 
1.3 % for a period of 20 days. 

In Figure 3 we compare the average radii from 1000 simulated data sets. We used a straight 
BW method and also the Fourier method. In the straight method we use a multidimentional 
minimization routine and a straight average of all 1000 radii. The mean radii of this method 
turn out to be significantly too large due to some outliers in the long tail of large radii. 
The Fourier method examines the errors and therefore averages the inverse radii for large 
photometric errors. The results of the Fourier method are very close to the true radius 
even for the largest photometric errors. We conclude that the straight BW radii are biased 
towards larger radii, especially for large photometric errors and short period stars. The 
Fourier method gives the same distribution of radii, but because the errors are understood, 
the mean radius is unbiased. 

The bias towards larger radii due to photometric errors in the straight BW method is depen- 
dent on A<pj and thus period dependent. We expect that the radii of short period Cepheids 
in a large sample will on average be overestimated, leading to a larger offset and smaller 
slope in the PR-relation. One way to reduce this systematic bias is to use an infrared (IR) 
bandpass (Welch 1994, Laney and Stobie 1995) for the light curve. A(f)j is larger if the 
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Figure 3: Mean radii for % velocity error, given photometric error, period and number of points. The 
straight BW radii are calculated using a multidimensional minimization routine. Shown is the average over 
1000 simulations. The error bars are the error in the mean. Note the change in scale from the straight BW 
radii to the Fourier radii. 

magnitude is measured in the IR, because the IR magnitude of the Cepheid is more sensitive 
to the radius than to the temperature. We thus expect that the slope of the PR relation 
is steeper, if IR magnitudes are used. This is confirmed by the results of Laney & Stobie 
(1995). 

3.2 Biases due to Truncation of the Fourier Series 

The small bias in the Fourier radii for no measurement errors in Figure 3 is due to truncation 
of the Fourier series. This bias increases if there are fewer points available. If there are only 
very few points then the maximum and the minimum of the velocity curve are underesti- 
mated. Thus the Fourier amplitudes become smaller. For our particular choice of velocity 
and light curves the Fourier amplitudes of the velocity curve are more underestimated than 
the amplitudes of the light curve. Therefore the average radius is estimated to be 0.3 % 
smaller than the true radius. This effect is very small compared to the measurement errors. 
But it suggests that the Fourier method is biased for data sets with a very small number of 
phase points. 

3.3 Biases due to Linearization 

For periods shorter than 20 days the effects of linearizing the logarithm are small. They are 
easily hidden in the measurement errors. In order to demonstrate that our scheme to correct 
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for the effects linearization works we simulated data sets with periods of 50 and 100 days. 
The radii for those periods are 159i? Q and 229i? Q and the ratios of At/Rq are 0.17 and 0.23 
respectively. 

We simulated 1000 data sets for each period with both 20 and 40 phase points. There 
were no deviations from the true radius inconsistent with the uncertainty of the mean. The 
correction factor was typically 2.2 % for a period of 50 days and 3.5 % for a period of 100 
days. The errors in the mean were much smaller than 1 % . We therefore conclude that our 
scheme to eliminate the effects of linerarization is successful. 

4 The Fourier Method Applied to Real Data 

To demonstrate our new method we applied it to a set of real Cepheid data. We chose the 
dataset of Bersier et al. (1994) because it provides measurement errors for each measurement. 
We measured the radii of the 19 Cepheids which had both photometry and velocities. 

We used all photometric points except for those of quality zero (see Bersier et al. 1994). We 
derived the radii using the V magnitudes, the B V color index and the radial velocities. The 
photometry is in the Geneva system, which is slightly different from the Johnson system, 
but has been used for BW solutions before (Burki & Meylan 1986). We assigned an error 
to each photometric point by the formula 0.02 / q where q is the quality of the photometry 
ranging from 1 to 4. We used the same relation for the quality of the color index, p, to get 
the errors of the color index. These formulae give errors for which the reduced x 2 of the fit 
is close to 1 for most objects. We used the velocity errors given by Bersier et al. (1994). 
For 3 of the stars, DX Gem, V465 Mon and W Sgr, we subtracted the velocity of the binary 
orbit before doing a Fourier fit. 

The number of Fourier orders we fitted for is again given by the largest phase gap. Comparing 
our Fourier coefficients to those derived by Bersier et al. (1994) we found that they are in 
good agreement. 

In Table 1 we show the number of data points per star and the quality of the Fourier fit. 
The Fourier method might not be adequate for any of those stars for which the reduced x 2 
of any of the three Fourier fits is large. The method might also not work for those stars 
whose data is only sufficient to fit 3 or less Fourier orders. Often the Cepheid curves cannot 
be well represented with 3 or less Fourier orders. The fits for SV Vul are particularly bad 
because the period is variable. The last column gives the mean residuals from a linear color 
index-surface brightness relation. Typical photometric errors are about 1 % . Thus most of 
the stars with good data seem to be consistent with a linear color index-surface brightness 
relation. 

The radii we found are shown in Table 2 and Figure 4. We used a constant p-factor of 
1.36 (Burki et al. 1982). The radius is proportional to the p-factor (eq. (2.8)), thus the 



13 



Table 1: Data and quality. n p hot and n ve i are the respective number of photometric and velocity data 
points. f ve i and f v hot are the number of Fourier orders fitted for. xlei' Xtemp and Xmag are the respective 
reduced x 2 values of the Fourier fits. The last column gives the mean deviation from a linear color index- 
surface brightness relation in magnitudes. 
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radii given in Table 2 can easily be converted to a different constant p-factor. The errors 
in radius are dominated by the uncertainty in the phase difference between color index and 
magnitude, A^j. 

4.1 Comparison of Radii 

Only some of the radii in Table 2 are well determined. Some data sets are not good enough 
to measure a believable radius with the Fourier method. DT Cyg for example has a Acfii 
that seems too small for its period. We fitted its light and color index curves with 4 Fourier 
orders instead of 3 and found that the resulting reduced x 2 of the fit is almost unchanged. 
However, A^i changes from 0.020 ± 0.007 to 0.034 ± 0.010 and the radius goes down from 
66.9 ~^28 2 8 8 R® t° 29.3 1iq2-R©. This shows that there are not enough photometric data 
points for DT Cyg to determine the radius with the Fourier method. 
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Figure 4: Period - Radius relation. 
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In Table 3 we compare the radii we have confidence in to previous radius determinations. 
The comparison shows that there are large differences in the various radii. The most overlap 
in data is between our radii and the radii by Ripepi et al. (1996), because both use the same 
velocity data. 

It is interesting to note that the PR-relation of Ripepi et al. (1996) is much shallower than 
the PR-relation of Gieren et al. (1989) or Laney & Stobie (1995). We also find a very 
shallow relation. The data Ripepi et al. (1996) and we have in common are the velocity 
data of Bersier et al. (1994). We suspect that this set of velocities requires the usage of a 
different p-factor to be consistent with the velocities used by other groups. 

As an example we examined the discrepancy between our result for T Vul and the result of 
Gieren, Barnes & Moffett (1989). We chose T Vul, because it shows the largest discrepancies. 
In order to exclude the effects of different BW realizations we applied our method to their 
data set. We used their B-V data and found A^i = 0.042 ± 0.003, in good agreement with 
our value. Thus the photometry is in agreement between the two data sets. The velocity 
curve of Gieren et al. (1989) is very noisy, leading to a first order Fourier amplitude of only 
11.7 km/sec compared to 14.4 km/sec for the data set of Bersier et al. (1994). This is the 
reason why we find a smaller radius of 39 -T^ 1 ^ Rq using the data set of Gieren et al. (1989). 
The discrepancy seems mostly due to noise in the velocity curve used by Gieren et al. (1989). 
The errors of the radii do overlap, thus they are not inconsistent. The 90 % confidence error 
bars on the Gieren et al. (1989) data is about 25 % . This partly explaines the large scatter 
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Table 2: Measurements of the mean radius of Cepheids, using a p factor of 1.36. A</>i is the phase difference 
between color index and magnitude in the first Fourier order. R is the corrected mean radius derived using 
eqs. (2.8). The errors in R are the 90 % confidence interval. The other errors are 1 a errors. c\ is the slope 
of the linear color index-surface brightness relation, obtained using eqs. (2.8). 
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in previous PR relations and the poor agreement among different measurements. 

A comparison with the maximum likelihood radii of Laney & Stobie (1995) is difficult because 
we do not have access to the data set used. It seems though like there is fairly good agreement 
after allowing for 10 % error in the radii of Laney & Stobie (1995). 

5 Conclusions 

We have demonstrated a new method of determining the radii of pulsating stars and their 
errors due to measurement uncertainties. By applying of our method to the dataset by Bersier 
et al. (1994) and by simulations we found that the uncertainty in each individual radius 
measurement using optical magnitudes is dominated by the error in the phase difference 
between color index and magnitude. This phase difference is smaller for shorter period 
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Table 3: Cepheid radii from recent BW solutions. GBM are the radii from Gieren, Barnes & Moffet (1989), 
RBMR the radii from Ripepi et al. (1996), LS V are the V, B-V radii from Laney & Stobie (1995) and LS 
K are the K, J-K radii of Laney & Stobie (1995). 
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Cepheids and leads to larger relative errors. This error is not caused by our Fourier method 
but affects any BW realization in the same way. Using the Fourier method we were able to 
identify this main source of measurement error. 

It is important to determine the radius error for each individual star. Often the errors are 
very large, explaining why different radius measurements of the same star arrive at very 
different results. The PR-relation changes if the errors are ignored, because the individual 
errors are often asymmetric and very different in size. 

If the measurement errors are known, the radius determination can be repeated on a large 
set of simulated data. The simulations will give the uncertainty and a possible bias due to 
the quality of the data and the method used. The radius determined from the data can then 
be corrected for this bias. 

It should be possible to eliminate all biases due to the particular BW method used with 
simulated data sets. Systematic biases due to the quality of the data can also be eliminated 
with simulations. A problem that still needs to be solved before reliable BW radii can be 
obtained is the value and the phase dependence of the p-factor. This will be addressed in a 
subsequent paper in this series. 
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